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The Weakly Interacting Massive Particle (WIMP) is one of the main candidates 
for the relic dark matter (DM). In the effective low-energy minimal supersymmetric 
standard model (effMSSM) the neutralino-nucleon spin and scalar cross sections in 
the low-mass regime were calculated. The calculated cross sections are compared 
<^ ■ with almost all experimental currently available exclusion curves for spin-dependent 

. WIMP-proton and WIMP-neutron cross sections. It is demonstrated that in general 

1 , about two-orders-of-magnitude improvement of the current DM experiment sensi- 

tivities is needed to reach the (effMSSM) SUSY predictions. At the current level 

; of accuracy it looks reasonable to safely neglect sub-dominant spin WIMP-nucleon 

contributions analyzing the data from spin-non-zero targets. To avoid misleading 

■ discrepancies between data and SUSY calculations it is, however, preferable to use 

a mixed spin-scalar coupling approach. This approach is applied to estimate future 
prospects of experiments with the odd-neutron high-spin isotope Ge. It is noticed 
that the DAMA evidence favors the light Higgs sector in the effMSSM, a high event 
rate in a 73 Ge detector and relatively high upgoing muon fluxes from relic neutralino 
annihilations in the Earth and the Sun. 

D 

' I. INTRODUCTION 



X 



Nowadays the main efforts and expectations in the direct dark matter searches are con- 
centrated in the field of so-called spin-independent (or scalar) interaction of a dark matter 
Weakly Interacting Massive Particle (WIMP) with a target nucleus. The lightest supersym- 
metric (SUSY) particle (LSP) neutralino is assumed to be the best WIMP dark matter (DM) 
candidate. It is believed that for heavy enough nuclei this spin-independent (SI) interaction of 
DM particles with nuclei usually gives the dominant contribution to the expected event rate of 
its detection. The reason is the strong (proportional to the squared mass of the target nucleus) 
enhancement of SI WIMP-nucleus interaction. The results currently obtained in the field are 
usually presented in the form of exclusion curves due to non-observation of the WIMPs. For 
a fixed mass of the WIMP the cross sections of SI elastic WIMP-nucleon interaction located 
above these curves are excluded. 



Only the DAMA collaboration claims observation of first evidence for the dark matter signal, 
due to registration of the annual modulation effect [1-3]. The DAMA results are shown in the 
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FIG. 1: WIMP-nucleon cross section limits in pb for scalar (spin-independent) interactions as a 
function of the WIMP mass in GeV. Shown are contour lines for some of the present experimental 
limits (solid lines) and some of projected experiments (dashed lines). The closed DAMA/Nal contour 
corresponds to a complete neglection of spin-dependent WIMP-nucleon interaction (<t sd = 0), while 
the open contour is obtained with the assumption that cr gD = 0.08 pb [2]. Our theoretical expectations 
are shown by scatter plots for a relic neutralino density 0.1 < ^ x ^o < 0.3 (g reen boxes) and to WMAP 
relic density 0.094 < fi^o < 0.129 (black triangles). Similar estimations one can find for example in 
[4-6]. 

middle of Fig. 1 as two contours together with some set of other exclusion curves already 
obtained (solid lines) and expected in the future (dashed lines). Aimed since more than one 
decade at the DM particle direct detection, the DAMA experiment (DAMA/Nal) with 100 
kg of highly radio-pure Nal(Tl) scintillator detectors successfully operated till July 2002 at 
the Gran Sasso National Laboratory of the I. N.F.N. On the basis of the results obtained over 
7 annual cycles (107731 kg-day total exposure) the presence of a WIMP model-independent 
annual modulation signature was demonstrated and the WIMP presence in the galactic halo is 
strongly supported at 6.3 a CL. [2]. The main result of the DAMA observation of the annual 




3 



modulation signature is the low-mass region of the WIMPs (40 < m x < 150 GeV), provided 
these WIMPs are cold dark matter particles. No other experiment at present has the sensitivity 
to look for this modulation effect. 

It is obvious that such a serious claim should be verified at least by one other completely 
independent experiment. To confirm this DAM A result one should perform a new experiment 
which would have (in reasonable time) the same or better sensitivity to the annual modula- 
tion signal (and also it would be better to locate this new setup in another low-background 
underground laboratory). This mission, in particular, could be executed by new-generation 
experiments with large enough mass of germanium HP detectors both with spin ( 73 Ge) and 
spin-less (natural Ge). Due to kinematic reasons (Mparget ~ -^wimp) these germanium isotopes 
with their masses being almost equal to the mass of the DAMA WIMP (about 70 GeV) have 
the best efficiency for such WIMP detection. A new setup with "naked" Ge detectors in liquid 
nitrogen (GENIUS-TF) is already installed and works over months under the low-background 
conditions of the Gran Sasso Laboratory [7]. The GENIUS-TF experiment is planned to be 
sensitive to the annual modulation signal with data taking over about 5 years with a large 
enough mass of the Ge detectors [8]. 

In this paper we start from the final results of the DAMA collaboration based on the 7-year- 
long measurements of the annual modulation [2] and consider their possible consequences for 
dark matter search with high-spin 73 Ge detectors like HDMS [9]. We also briefly consider some 
aspects of the spin-dependent (or axial-vector) interaction of the DM WIMPs with nuclei. 
There are at least three reasons to think that this SD interaction could be very important. 
First, contrary to the only one constraint for SUSY models available from the scalar WIMP- 
nucleus interaction, the spin WIMP-nucleus interaction supplies us with two such constraints 
(see for example [10] and formulas below). Second, one can notice [4, 11] that even with a very 
sensitive DM detector (say, with a sensitivity of 10~ 5 events/day/kg) which is sensitive only to 
the WIMP-nucleus scalar interaction (with spin-less target nuclei) one can, in principle, miss 
a DM signal. To safely avoid such a situation one should have a spin-sensitive DM detector, 
i.e. a detector with spin-non-zero target nuclei. Finally, there is a complicated nuclear spin 
structure, which for example, possesses the so-called long g-tail form-factor behavior. Therefore 
for heavy mass target nuclei and heavy WIMP masses the SD efficiency to detect a DM signal 
is much higher than the SI efficiency [12]. 
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II. APPROACH TO OUR CALCULATIONS 
A. Cross sections and event rate 

A dark matter event is elastic scattering of a relic neutralino x from a target nucleus A 
producing a nuclear recoil En which can be detected by a suitable detector. The differential 
event rate in respect to the recoil energy is the subject of experimental measurements. The 
rate depends on the distribution of the relic neutralinos in the solar vicinity f(v) and the cross 
section of neutralino-nucleus elastic scattering [2, 13-20]. The differential event rate per unit 
mass of the target material has the form 

dR AT p Y / >1 ' max , ,/ N da . 9 , 

-Nt—\ dvf(v)v—(v,q 2 ). (1) 



dE R m x J Vmin dq 

The nuclear recoil energy E R = q 2 /(2MA) is typically about lCT 6 m x and N T = M jA is the 
number density of target nuclei, where M is the Avogadro number and A is the atomic mass of 
the nuclei with mass Ma- The neutralino-nucleus elastic scattering cross section for spin-non- 
zero (J 7^ 0) nuclei contains SI and SD terms [12, 21, 22]: 

d ° A ( „„2, EI-MI 2 SjM) SiM 2 ) 



dq 2K ^' tvv 2 (2J+1) v 2 (2J+l) v 2 (2J+l) 

_ °"sd(0) p 2 ( 2x , gjilW p2 / 2\ / 9 \ 

The normalized (-^sdsi(O) — 1) non-zero-momentum-transfer nuclear form-factors 

FsDSli " } - (3) 

are defined via nuclear structure functions [12, 21, 22] 

S&q) = E |(J||C L (g)||J)r^|(J||C (g)||J)| 2 , (4) 

L even 

S&(q) = E m\\T£ l5 (q)\\N)\ 2 + \(N\\£l(q)\\N)f). (5) 

Lodd 

The transverse electric T el5 (q) and longitudinal C 5 (q) multipole projections of the axial vector 
current operator, and the scalar function Cl((?) are given in the form 

^ = ^TT?(T + 5(^)[^TTM i , i+1 ( 5 f i ) + VlM i , i _ 1 ( 9 r ( ) 
C h{q) = E c oJL(qri)Y L (ri), C (q) = E c jo(^)>o(^), 

i, nucleons « 
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where a Q1 = a n ±a p and M LjL >(qfi) = jL'('?^)[>L'(^)o ! j] L [12, 21, 22]. The nuclear SD and SI 
cross sections at q = in (2) can be presented as follows 



. 4/4Ssi(0)_/4 2 „ 



«&«>) = M^M=M(£±i){ ap( s-) + a „(s^} 2 



^ 3 J 
Atj4 J+l 



(7) 

^(0)(S p A ) + sign(a p a n )^ib70)(S^)}' (8) 
^(0){(S p A )cos^+(S T f)sin^} 2 . (9) 



Here ha = — rr~ is the reduced x~A mass, and a v = u n is assumed. Following Bernabei et 

m x + Ma 

al. [2, 23] the effective spin WIMP-nucleon cross section Osr>(0) and the coupling mixing angle 
6 were introduced 



^d(O) = + tan^ = ^; (10) 



7T 3 



^ D = a^-cos 2 £, a« D = all ■ sin 2 6. (11) 

The zero-momentum-transfer proton and neutron SI and SD cross sections 

<(0)=4^4 c^c^ECJ^); (12) 

" <? 

ct§?(0) = 12^<£ B a p = Y,AAf, a n = Y,A^ q n) (13) 

71 q q 

depend on the effective neutralino-quark scalar C q and axial-vector A q couplings from the 
effective Lagrangian 

£efr = ^2(Aq- xi»i5X ■ qYi$q + C q ■ xx ■ qq) + ••• (14) 

q 

and on the spin (A q P^) and mass (fjf' 71 ^) structure of nucleons. The parameters a p(n ) in (13) 
can be considered as effective WIMP-proton (neutron) couplings. The factors A q p ^ in (13) 
parameterize the quark spin content of the nucleon and are defined by the relation 2A q n ' p ^s^ = 
(p, s|'0q7' i 75'0 ? |p, s )(p,n)- A global QCD analysis for the g 1 structure functions [24] including 
0(al) corrections supplied us with the values [25] 

A<?) = Aj° = 0.78 ± 0.02, A? = A*"' = -0.48 ± 0.02, A® = Af' = -0.15 ± 0.02. (15) 

The nuclear spin (proton, neutron) operator is defined as follows 

A 

S Pl „ = EvW. (16) 
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where % runs over all nucleons. Further the convention is used that all angular momentum 
operators are evaluated in their z-projection in the maximal Mj state, e.g. 

(S) = (N\S\N) ee (J, Mj = J\S z \J,Mj = J). (17) 

Therefore (S p („)) is the spin of the proton (neutron) averaged over all nucleons in the nucleus A. 
The cross sections at zero momentum transfer show strong dependence on the nuclear structure 
of the ground state [26-28]. 

The relic neutralinos in the halo of our Galaxy have a mean velocity of (v) ~ 300 km/s = 
10~ 3 c. When the product q max R <C 1, where R is the nuclear radius and g max = 2[iav is 
the maximum momentum transfer in the x~A scattering, the matrix element for the SD \-A 
scattering reduces to a very simple form {zero momentum transfer limit) [26, 27]: 

M = C{N\a p S p + a n S n \N) ■ s x = CA(N\ J\N) ■ s x . (18) 

Here s x is the spin of the neutralino, and 

(N\a p S p + a n S n \N) _ (N\(a p S p + a n S n ) ■ 3\N) 

(N\J\N) J(J+1) ' 1 ] 

It is seen that the x couples to the spin carried by the protons and the neutrons. The nor- 
malization C involves the coupling constants, masses of the exchanged bosons and various LSP 
mixing parameters that have no effect upon the nuclear matrix element [29]. In the q = limit 
the spin structure function (5) reduces to 

2 7-1-1 

Si D (0) = ^—A 2 J(J+l). (20) 

71 

The first model to estimate the spin content in the nucleus for the dark matter search was 
the independent single-particle shell model (ISPSM) used originally by Goodman and Witten 
[30] and later in [15, 31, 32]. Here the ground state value of the nuclear total spin J can be 
described by that of one extra nucleon interacting with the effective potential of the nuclear 
core. There are nuclear structure calculations (including non-zero-momentum approximation) 
for spin-dependent neutralino interaction with helium 3 He [33]; fluorine 19 F [28, 33, 34]; sodium 
23 Na [27, 28, 33, 34]; aluminum 27 M [26]; silicon 29 Si [22, 28, 34]; chlorine 35 C1 [22]; potassium 
39 K [26]; germanium 73 Ge [22, 35]; niobium 93 Nd [36]; iodine 127 I [27]; xenon 129 Xe [27] and 
131 Xe [12, 27, 37]; tellurium 123 Te [37] and 125 Te [27]; lead 208 Pb [33, 38]. The zero-momentum 
case is also investigated for Cd, Cs, Ba and La in [37, 39, 40]. 

There are several approaches to more accurate calculations of the nuclear structure effects 
relevant to the dark matter detection. The list of the models includes the Odd Group Model 
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(OGM) of Engel and Vogel [41] and their extended OGM (EOGM) [21, 41]; Interacting Bo- 
son Fermion Model (IBFM) of Iachello, Krauss, and Maino [40]; Theory of Finite Fermi Sys- 
tems (TFFS) of Nikolaev and Klapdor-Kleingrothaus [42]; Quasi Tamm-Dancoff Approximation 
(QTDA) of Engel [12]; different shell model treatments (SM) by Pacheco and Strottman [39]; 
by Engel, Pittel, Ormand and Vogel [36] and Engel, Ressell, Towner and Ormand, [26], by 
Ressell et al. [22] and Ressell and Dean [27]; by Kosmas, Vergados et al. [28, 33, 38]; the 
so-called "hybrid" model of Dimitrov, Engel and Pittel [35] and perturbation theory based on 
calculations of Engel et al. [26]. 

The direct detection rate (1) in a nucleus A integrated over the recoil energy interval from 
threshold energy, e, till maximal energy, e, is a sum of SD and SI contributions: 



R(e,e) = a(e,e,mj<7jg 1 + /3(e,e,77i x )ajg£; (21) 
a(e,e,m x ) = N T ^^A 2 A sl (e,e), 



P(e,e,m x )=N T ^^^^((S^cos9+(S^sm9) 2 A SD (e,e)- 

A S i,s D (e, e) = H jT dE R Fl SB (E R )I(E R ). (22) 

To estimate the event rate (21) one needs to know a number of quite uncertain astrophysical 
and nuclear structure parameters as well as the precise characteristics of the experimental setup 
(see, for example, the discussions in [2, 43]). 

B. Effective low-energy MSSM 

To obtain as much as general predictions it appeared more convenient to work within a 
phenomenological SUSY model whose parameters are defined directly at the electroweak scale, 
relaxing completely constraints following from any unification assumption as for example in 
[44-49], and which is called an effective scheme of MSSM (effMSSM) in [50], and later by 
some people low-energy effective supersymmetric theory (LEEST) in [5, 6]. In our previous 
calculations in effMSSM [4, 11, 18-20, 44, 51-54] we have adopted some effective scheme (with 
non-universal scalar masses and with non-universal gaugino soft masses) which lead to large 
values for direct detection rates of DM neutralinos. 

Our MSSM parameter space is determined by the entries of the mass matrices of neutralinos, 
charginos, Higgs bosons, sleptons and squarks. The relevant definitions one can find in [44] . The 
list of free parameters includes: tan (3 is the ratio of neutral Higgs boson vacuum expectation 
values, n is the bilinear Higgs parameter of the superpotential, Mi j2 are soft gaugino masses, Ma 
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is the CP-odd Higgs mass, m~, m~, m~ (rri~, m~) are squared squark (slepton) mass parameters 

CJ U U L, ill 

for the 1st and 2nd generation, m~ , m~, m~ (m~ , m~) are squared squark (slepton) mass 
parameters for 3rd generation and A t , A^, A T are soft trilinear couplings for the 3rd generation. 
The third gaugino mass parameter M3 defines the mass of the gluino in the model and is 
determined by means of the GUT assumption M2 = 0.3 M3. 

Contrary to our previous considerations [4, 11, 18-20, 44, 51-54] and aiming at exploration 
of the MSSM parameter space in the DAMA-inspired domain of the lower masses of the LSP 
(m x < 200 GeV), we narrowed in the present work the intervals of the randomly scanned 
parameter space to the following: 

-200 GeV < Mi < 200 GeV, -1 TeV < M 2 , fi < 1 TeV, -2 TeV < A t < 2 TeV, 

10 < tan/? < 50, 50 GeV <M A < 500 GeV, (23) 
10 GeV 2 < m 2 ~~ ,m 2 ~~ < 10 6 GeV 2 . 

As previously we assume that squark masses are basically degenerate. Bounds on flavor- 
changing neutral currents imply that squarks with equal gauge quantum numbers must be 
close in mass [55-57]. With the possible exception of third generation squarks the assumed 
degeneracy holds almost model-independently [55]. Therefore for other sfermion mass param- 
eters as before in [4, 11, 44, 51-54] we used the relations m~ = m~ = m~, m 2 ~ = m~, 

U D (eg Cj Li 

m% = m 2 ~ = m~ , m~ = m~ . The parameters At, and A T are fixed to be zero. 

We have included the current experimental upper limits on sparticle and Higgs masses from 
the Particle Data Group [58]. For example, we use as previously the following lower bounds 
for the SUSY particles: M^± > 100 GeV for charginos, M^ 23 > 45,76,127 GeV for non- 
LSP neutralinos, respectively; M„ > 43 GeV for sneutrinos, M- eR > 70 GeV for selectrons, 
Mq > 210 GeV for squarks, M h > 85 GeV for light top-squark, M H o > 100 GeV for neutral 
Higgs bosons, M H + > 70 GeV for the charged Higgs boson. Also the limits on the rare b — > 57 
decay [59, 60] following [61-64] have been imposed. 

For each point in the MSSM parameter space (MSSM model) we have evaluated the relic 
density of the light neutralinos Q x hl with our code [52-54] based on [65], taking into account 
all coannihilation channels with two-body final states that can occur between neutralinos, 
charginos, sleptons, stops and sbottoms, as long as their masses are rrii < 2m x . We assume 
as before 0.1 < VL x h 2 < 0.3 for the cosmologically interesting region and we also consider the 
WMAP reduction of the region to 0.094 < Vt x h 2 < 0.129 [66, 67] and a possibility the LSP to 
be not a unique DM candidate with much smaller relic density 0.002 < Vth 2 < 0.1. 
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III. RESULTS AND DISCUSSIONS 
A. Cross sections in the effMSSM for m x < 200 GeV 

The results of our evaluations of the zero-momentum-transfer proton and neutron SI (12) 
and SD (13) cross sections in the effMSSM approach within the DAMA-inspired parameter 
space of (23) are shown as scatter plots in Figs. 2-4. 




50 li'i'.i 150 50 1flf> 150 



FIG. 2: Cross sections of the spin-dependent (spin) and the spin-independent (scalar) interactions 
of WIMPs with the proton and the neutron. Filled green circles correspond to the relic neutralino 
density < Ov/iq < 1, red squares correspond to the sub-dominant relic neutralino contribution 
0.002 < O x /tg < 0.1 and black triangles correspond to the relic neutralino density 0.1 < Q-yho < 0.3 
(left panel) and to the WMAP relic density 0.094 < Q x hl < 0.129 (right panels). 

Scatter plots with individual cross sections of SD and SI interactions of WIMPs with proton 
and neutron are given in Fig. 2 as functions of the LSP mass. In the figure filled green circles 
correspond to cross sections calculated when the neutralino relic density should just not over- 
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close the Universe (0.0 < £l x hl < 1.0). Filled red squares show the same cross sections when 
one assumes the relic neutralinos to be not the only DM particles and give only sub-dominant 
contribution to the relic density 0.002 < VL x h^ < 0.1. In the left panel of Fig. 2 these cross 
sections are shown with the black triangles corresponding to the case when the relic neutralino 
density is in the bounds previously associated with a so-called flat and accelerating Universe 
0.1 < Q x hl < 0.3. The black triangles in the right panel in Fig. 2 correspond to imposing 
the new WMAP [66, 67] constraint on matter relic density 0.094 < £l x hl < 0.129. Despite a 
visible reduction of the allowed domain for the relic density due to the WMAP result the upper 
bounds for the spin-dependent and the spin-independent WIMP-nucleon cross section are not 
significantly affected. From the comparison of circle and square distributions, as expected, 
follows that the largest cross section values correspond to smallest values of the Q x , especially 
for smaller LSP masses. It is seen that the LSP as a sub-dominant DM particle favors the large 
SD and SI cross sections. Furthermore the maximal SD and SI cross sections in Fig. 2 (green 
circles) come for very small relic density values 0.0 < Q x hl < 0.002. 

One can see also that in our efiMSSM with parameters from (23) the lower bound value in 
the relic density constraints (as for example, 0.094 in the case of WMAP) restricts from below 
the allowed masses of the LSP in accordance with previous considerations [18, 68]. 

The spin-dependent and spin-independent WIMP-proton cross sections as functions of input 
MSSM parameters /x, M A) tan/5 and m~ are shown in Figs. 3 and 4. There is no noticeable 
dependence of these scatter plots on the other free parameters from our set (23). From these 
figures one can see the similarity of the scatter plots for spin-dependent and and scalar cross 
sections as functions of /i and m~. Decrease of both lower bounds of the cross sections with m~ 
occur due to increase of masses of squarks, which enter the s-channel intermediate states. Both 
spin-dependent and spin-independent cross sections increase when decreases, in agreement 
with literature [45, 69, 70] and our previous calculations [4, 44]. The increase of the scalar 
cross sections generally is connected with an increase of the Higgsino admixture of the LSP 
and increase of Higgsino-gaugino interference which enters this cross section [69-71]. The 
reason of the Higgsino growth can be non- universality of scalar soft masses [70], variation of 
intermediate unification scale [71] , or focus point regime of the supersymmetry [69] . There is no 
any visible sensitivity of the SD cross sections to tan (3 and Ma (Higgs bosons do not contribute) 
but the SI cross section possesses remarkable dependence on these parameters. The SI cross 
sections rather quickly drop with growth of the CP-odd Higgs mass M A and increase with tan (3 
[20, 44, 50, 70-73]. The different tan (3- and M^-dependence of the SD and SI cross section 
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FIG. 3: Cross sections of of WIMP-proton spin-independent interactions as function of input param- 
eters fi, Ma, tan (3 and m~ with the same notations as in Fig. 2 for used constraints on the neutralino 

relic density Q x h 2 . 

as well as the general about-2-order-of-magnitude excess of the spin-dependent cross sections 
over spin-independent cross sections may be important for observations [4, 11, 51, 74, 75]. It is 
interesting to note, that maximal values for the LSP-proton SD cross section one can obtain in 
the pure Higgsino case (when only Z-exchange contributes) at a level of 5 ■ 10~ 2 pb. This value 
is almost reached by points from our scatter plots in Figs. 2 and 4. 

B. Constraints on WIMP-nucleon spin interactions 

For the spin-zero nuclear target the experimentally measured event rate (1) of direct DM par- 
ticle detection, via formula (2) is connected with the zero-momentum WIMP-proton(neutron) 
cross section (6). The zero momentum scalar WIMP-proton(neutron) cross section crgj(O) can 
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FIG. 4: Cross sections of of WIMP-proton spin-dependent interactions as function of input parameters 
fj,, Ma, tan [5 and m~ with the same notations as in Fig. 2 for used constraints on the neutralino relic 

density fl x h 2 . 

be expressed through effective neutralino-quark couplings C q (14) by means of expression (12). 
These couplings C q (as well as A q ) can be directly connected with the fundamental parameters 
of a SUSY model such as tan/3, M 12 , /i, masses of sfermions and Higgs bosons, etc. Therefore 
experimental limitations on the spin-independent neutralino-nucleon cross section supply us 
with a constraint on the fundamental parameters of an underlying SUSY model. 

In the case of the spin-dependent WIMP-nucleus interaction from a measured differential 
rate (1) one first extracts a limitation for 0"^>(O), and therefore has in principle two constraints 
[10] for the neutralino-proton a p and neutralino-neutron a n effective spin couplings, as follows 
from relation (7). From (7) one can also see that contrary to the spin-independent case (6) 
there is no, in general, factorization of the nuclear structure for cr^ D (0). Both proton (S^ ) and 
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neutron (S^) spin contributions simultaneously enter into formula (7) for the SD WIMP-nucleus 
cross section crg D (0). 

In the earlier considerations based on the OGM [21, 41] one assumed that the nuclear spin 
is carried by the "odd" unpaired group of protons or neutrons and only one of either (S„) or 
(Sp) is non-zero (the same is true in the ISPSM [15, 30-32]). In this case all possible target 
nuclei can naturally be classified into n-odd and p-odd groups. 

Following this classification the current experimental situation in the form of the exclusion 
curves for the spin-dependent WIMP-proton cross section is given in Fig. 5. The data are 
taken from experiments BRS, (Nal, 1992) [76, 77], BPRS (CaF 2 , 1993) [78], EDELWEISS 
(sapphire, 1996) [79], DAMA (NAI, 1996) [80], DAMA (CaF 2 , 1999) [81, 82], UKDMS (Nal, 
1996) [83-86], ELEGANT (CaF 2 , 1998) [87], ELEGANT (Nal, 1999) [88, 89], Tokio (LiF, 1999, 
2002) [90-94], SIMPLE (C 2 C1F 5 , 2001) [95], CRESST (A1 2 3 , 2002) [96], PICASSO (C n F m , 
2002) [97], ANAIS (Nal, 2002) [98] and NAIAD (Nal, 2003) [99]. Although the DAMA/NaI-7 
(2003) contours [2] are obtained on the basis of the positive and model-independent signature 
of the annual signal modulation (closed contour) as well as in the mixed coupling framework 
(open contour) [23] the contours for the WIMP-proton SD interaction (dominating in 127 I) are 
also presented in the figure (we will discuss the situation later). 

The current experimental situation for the spin-dependent WIMP-neutron cross sections 
is given in Fig. 6. The data are taken from the first experiments with natural Ge (1988, 1991) 
[100, 101], xenon (DAMA/Xe-0,2) [102-104], sodium iodide (NAIAD) [99], and from the HDMS 
experiment with a 73 Ge target [105]. Similar to Fig. 5 the DAMA/NaI-7 (2003) [2] contours 
for the WIMP-neutron SD interaction (sub-dominant in 127 I) are placed in the figure. In the 
future one can also expect some exclusion curves for the SD cross section, for example, from the 
CDMS [106] and EDELWEISS [107] experiments with natural germanium bolometric detectors. 

In Figs. 5 and 6 are also given scatter plots for SD proton and neutron cross sections which 
correspond to the results of our calculations shown in Figs. 2-4. From Figs. 5 and 6 one 
can, in general, conclude that an about two-orders-of-magnitude improvement of the current 
DM experimental sensitivities (in the form of these exclusion curves) is needed to reach the 
SUSY predictions for the of^ , provided the SUSY lightest neutralino is the best WIMP particle 
candidate. 

Here we note that the calculated scatter plots for <jg D (Fig. 5) are obtained without any 
assumption about o"5 D = 0, but the experimental exclusion curves for Og D traditionally were 
extracted from the data under the full neglection of the spin-neutron contribution, i.e. under 
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FIG. 5: Currently available exclusion curves for spin-dependent WIMP-proton cross sections (cgrj as 
a function of WIMP mass). The curves are obtained from [76-99]. DAMA/NaI-7a(f) contours for 
WIMP-proton SD interaction in 127 I are obtained on the basis of the positive and model independent 
signature of annual signal modulation in the framework of a mixed scalar-spin coupling approach 
[2, 23]. The scatter plots correspond to our calculations given in Figs. 2-4. The small triangle-like 
shaded area in the bottom is taken from [25]. Note that the closed DAMA contour is above the upper 
limit for <jg D 5 • 1(T 2 pb. 



the assumption crg D = 0. This one-spin-coupling dominance scheme (always used before a 
new approach was proposed in [74]) gave a bit too pessimistic exclusion curves, but allowed 
direct comparison of exclusion curves from different experiments. More stringent constraints 
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FIG. 6: Currently available exclusion curves for spin-dependent WIMP-neutron cross sections (crg D 
versus WIMP mass). The curves are taken from [99-105]. DAMA/NaI-7a(f) contours for WIMP- 
neutron SD interaction (sub-dominating in 127 1) are obtained by us from the relevant figures of [2, 23]. 
The scatter plots correspond to our calculations given in Figs. 2-4. Note that the NAIAD curve here 
corresponds to the sub-dominant for 127 I WIMP-neutron SD interaction. The curve was extracted 
from the nucleus 127 I (which has dominating WIMP-proton SD interaction) in the approach of [74]. It 
is much weaker in comparison with the relevant NAIAD curve for the WIMP-proton SD interaction 
in Fig. 5. 

on <7g D one obtains following [74] and [1-3] assuming both <7g D ^ and <7g D ^ although 
usually the contribution of the neutron spin is very small ((S^) <C (S^ )). Therefore the direct 
comparison of old-fashioned exclusion curves with new ones is misleading in general. The same 
conclusion concerns [2, 3] direct comparison of the SI exclusion curves (obtained without any 
SD contribution) with new SI exclusion curves (obtained with non-zero SD contribution) as 
well as with the results of the SUSY calculations (Fig. 1). 
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C. Mixed spin-scalar WIMP-nucleon interactions 

Further more accurate calculations of spin nuclear structure [12, 22, 26-28, 33, 35, 36, 38- 
40] demonstrate that contrary to the simplified odd-group approach both (S p ) and (S^) differ 
from zero, but nevertheless one of these spin quantities always dominates ((Sp 4 ) <C (S^), or 
(S£) < (S^)). // together with the dominance like (S^ (n) ) < (S^ (p) ) one would have the WIMP- 
proton and WIMP-neutron couplings of the same order of magnitude (not a„( p ) <^a p ^), the 
situation could look like that in the odd-group model and one could safely (at the current level 
of accuracy) neglect sub-dominant spin contribution in the data analysis. 

Nevertheless it was shown in [74] that in the general SUSY model one can meet right a 
case when a n ( p ) <C a p(n ) and proton and neutron spin contributions are strongly mixed. To 
separately constrain the SD proton and neutron contributions at least two new approaches 
appeared in the literature [23, 74]. As the authors of [74] claimed, their method has the 
advantage that the limits on individual WIMP-proton and WIMP-neutron SD cross sections 
for a given WIMP mass can be combined to give a model-independent limit on the properties of 
WIMP scattering from both protons and neutrons in the target nucleus. The method relies on 
the assumption that the WIMP-nuclear SD cross section can be presented in the form a^. (0) = 
V^U) 2 . where ohU and o SdU « ^ — es, not dne ctly colled 
with measurements. Furthermore, to extract, for example, a constraint on the sub-dominant 
WIMP-proton spin contribution one should assume the proton contribution dominance for a 
nucleus whose spin is almost completely determined by the neutron-odd group. From one side, 
this may look almost useless, especially because these sub-dominant constraints are always 
much weaker than the relevant constraints obtained directly with a proton-odd group target 
(one can compare, for example, the restrictive potential of the NAIAD exclusion curves in 
Figs. 5 and 6). From another side, the very large and very small ratios o v ja n ~ a p /a n obtained 
in [74] correspond to neutralinos which are extremely pure gauginos. In this case Z-boson 
exchange in SD interactions is absent and only sfermions give contributions to SD cross sections. 
Obviously this is a very particular case which is also currently not in agreement with the 
experiments. Following an analogy between neutrinos and neutralinos on could assume that 
neutralino couplings with the neutron and the proton should not to very different [108] and one 
could expect preferably |a n |/|o p | ~ 0(1). We have checked the assumption in our effMSSM 
approach for large LSP masses in [51, 52] and in this paper for relatively low LSP masses 
m x < 200 GeV. Figure 7 shows that for the ratio of a n to a p we have the bounds 
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FIG. 7: The scatter plots (circles, squares and triangles) give the ratio of the neutralino-neutron 
spin coupling a n to the neutralino-proton spin coupling a p in the effMSSM under the notations as in 
Figs. 2-4. The ratio is restricted to the range between 0.55 and 0.8. 



0.55 < 



< 0.8. 



(24) 



The scatter plots in Fig. 7 as previously (see Fig. 2) were obtained with the relic neutralino 
density 0.0 < ^l x h^ < 1.0 (green circles), with sub-dominant relic neutralino contribution 
0.002 < f2 x /ig < 0.1 (red squares) and with a WMAP-inspired relic neutralino density of 
0.094 < {l x h,Q < 0.129 (black triangles). Therefore in the model the couplings are almost the 
same and one can safely neglect, for example, the (S^ )-spin contribution in the analysis of the 
DM data for a nuclear target with (S£) <C (S^). 

Furthermore, when one compares in the same figure an exclusion curve for SD WIMP- 
proton coupling obtained without sub-dominant SD WIMP-neutron contribution and without 
SI contribution (all curves in Fig. 5 except the one for NAIAD [99] and one for Tokyo-LiF [94]), 
with a curve from the approach of [74], when the sub-dominant contribution is included (the 
NAIAD and Tokyo-LiF curves in Fig. 5), one "artificially" improves the sensitivity of the latter 
curves (NAIAD or Tokyo-LiF) in comparison with the former ones. To be consistent and for 
reliable comparisons, one should coherently recalculate all previous curves in the new manner. 
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This message is clearly also stressed in [2]. The same arguments are true for the last results 
of the SIMPLE experiment [109] and search for DM with NaF bolometers [110] where the 
SI contribution seems also completely ignored. Both above-mentioned results for fluorine will 
obviously be worse if (contrary to calculations of [39]) more reliable 19 F spin matrix elements 
(for example, from [28]) were used in their analysis. Although 19 F has the best properties 
for investigation of WIMP-nucleon spin-dependent interactions (see, for example [28]) it is not 
obvious that one should completely ignore spin-independent WIMP coupling with the fluorine. 

which follows from (6)-(9), it is 



A,p 



For example, in the relation a ~ a s ^ 



not a priori clear that <C At least for isotopes with atomic number A > 50 [13, 20] 

CT SD CT SD 

the neglection of the SI contribution would be a larger mistake than the neglection of the sub- 
dominant SD WIMP-neutron contribution, when the SD WIMP-proton interaction dominates. 
Therefore we would like to note that the "old" odd-group-based approach in analyzing the 
SD data from experiments with heavy enough targets (for example, germanium) is still quite 
suitable. Especially when it is not obvious that (both) spin couplings dominate over the scalar 
one. 

From measurements with 73 Ge one can extract, in principle, not only the dominant constraint 
for WIMP-nucleon coupling a n (or <7g D ) but also the constraint for the sub-dominant WIMP- 
proton coupling a p (or <7g D ) using the approach of [74]. Nevertheless, the latter constraint will 
be much weaker in comparison with the constraints from p-odd group nuclear targets, like 19 F 
or I. This fact is illustrated by the NAIAD (Nal, 2003) curve in Fig. 6, which corresponds to 
the sub-dominant WIMP-neutron spin contribution extracted from the p-odd nucleus I. 

Another approach of Bernabei et al. [23] looks in a more appropriate way for the mixed 
spin-scalar coupling data presentation, and is based on an introduction of the so-called effective 
SD nucleon cross section o"sd(0) (originally <tsd in [2, 23]) and coupling mixing angle 9 (10) 
instead of crg D (0) and <7g D (0). With these definitions the SD WIMP-proton and WIMP-neutron 
cross sections are given by relations (11). 

In Fig. 8 the WIMP-nucleon spin and scalar mixed couplings allowed by the annual modu- 
lation signature from the 100- kg DAMA/Nal experiment are shown inside the shaded regions. 
The regions from [2, 23] in the (£o"si, £°sd) space for 40 GeV< rn wlMP <H0 GeV cover spin- 
scalar mixing coupling for the proton (6 = case of [2, 23], left panel) and spin-scalar mixing 
coupling for the neutron (6 = 7r/2, right panel). From nuclear physics one has for the proton 
spin dominated 23 Na and 127 I < 0.1 and < 0.02^-0.23, respectively. For the 9 = due 
to the p-oddness of the I target, the DAMA WIMP-proton spin constraint is the most severe 
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FIG. 8: The DAMA-Nal region from the WIMP annual modulation signature in the (£o"si, £o"sd) 
space for 40 < rn wlMP < 110 GeV [2, 23]. Left panel corresponds to dominating (in 127 I) SD proton 
coupling only (9 = 0) and right panel corresponds to sub-dominating SD neutron coupling only {6 = 
7r/2). The scatter plots give correlations between Cgj and <r SD in the effMSSM (£ = 1 is assumed) for 
m x < 200 GeV under the same notations as in Figs. 2-4. In the right panel the DAMA liquid xenon 
exclusion curve from [23] is given (dashed line). 

one (see Fig. 5). 

In the right panel of Fig. 8 we present the exclusion curve (dashed line) for the WIMP- 
proton spin coupling from the proton-odd isotope 129 Xe obtained under the mixed coupling 
assumptions [23] from the DAMA-LiXe (1998) experiment [104, 111, 112]. For the DAMA Nal 
detector the 8 = tt/2 means no (S p ) contribution at all. Therefore, in this case DAMA gives 
the sub-dominant (S n ) contribution only, which could be compared further with the dominant 
(S n ) contribution in 73 Ge. 

The scatter plots in Fig. 8 give Ugj as a function of <7g D (left panel) and <Tg D (right panel) 
calculated in this work in the effMSSM with parameters from (23) under the same constraints 
on the relic neutralino density as in Figs. 2-4. Filled circles (green) correspond to relic neu- 
tralino density 0.0 < fi x h^ < 1.0, squares (red) correspond to sub-dominant relic neutralino 
contribution 0.002 < fi x /ig < 0.1 and triangles (black) correspond to WMAP density constraint 
0.094 < tt x h 2 < 0.129. 

The constraints on the SUSY parameter space in the mixed coupling framework in Fig. 8 
are, in general, much stronger in comparison with the traditional approach based on the one- 
coupling dominance (Figs. 1, 5 and 6). 

It follows from Fig. 8, that when the LSP is the sub-dominant DM particle (squares in the 
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figure) SD WIMP-proton and WIMP-neutron cross sections at a level of 3 5 • 10~ 3 pb are 
allowed, but the WMAP relic density constraint (triangles) together with the DAMA restric- 
tions leaves only < 3 • 10~ 5 pb without any visible reduction of allowed values for <7gj. In 
general, according to the DAMA restrictions, small SI cross sections are completely excluded, 
only o"si > 3 4- 5 ■ 10~ 7 pb are allowed. Concerning the SD cross section the situation is not 
clear, because for the allowed values of the SI contribution, the SD DAMA sensitivity did not 
yet reach the calculated upper bound for the SD LSP-proton cross section of 5 • 10~ 2 pb. 

D. The mixed couplings case for the high-spin 73 Ge 

Comparing the number of exclusion curves in Figs. 5 and 6 one can see that there are many 
measurements with p-odd nuclei and there is a lack of data for n-odd nuclei, i.e. for <jg D . 
Therefore measurements with n-odd nuclei are needed. From our point of view this lack of 
<7g D measurements can be filled with new data expected from the HDMS experiment with the 
high-spin isotope 73 Ge [9]. This isotope looks with a good accuracy like an almost pure n-odd 
group nucleus with (S n ) (S p ) (Table I). The variation of the (S p ) and (S n ) in the table 
reflects the level of inaccuracy and complexity of the current nuclear structure calculations. 

TABLE I: All available calculations in different nuclear models for the zero-momentum spin structure 
(and predicted magnetic moments /i) of the Ge nucleus. The experimental value of the magnetic 
moment given in the brackets is used as input in the calculations. 



73 Ge (G 9/2 ) 


<s P ) 


(Sn) 


\i (in ix N ) 


ISPSM, Ellis-Flores [32, 113] 





0.5 


-1.913 


OGM, Engel-Vogel [41] 





0.23 


(-0.879) cxp 


IBFM, Iachello at al. [40] and [22] 


-0.009 


0.469 


-1.785 


IBFM (quenched), Iachello at al. [40] and [22] 


-0.005 


0.245 


(-0.879) cxp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [42] 





0.34 




SM (small), Ressell at al. [22] 


0.005 


0.496 


-1.468 


SM (large), Ressell at al. [22] 


0.011 


0.468 


-1.239 


SM (large, quenched), Ressell at al. [22] 


0.009 


0.372 


(-0.879) cxp 


"Hybrid" SM, Dimitrov at al. [35] 


0.030 


0.378 


-0.920 



In the mixed spin-scalar coupling case for 73 Ge the direct detection rate integrated over 
recoil energy (21) from threshold energy, e, till maximal energy, e can be presented in the form 

R(e,e) = a(e,e,m x )a§ I + /3(e,e,m x )<7g D ; (25) 
a(e,e,m x ) = N T p^A 2 A sl (e,e), 
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P(e,e,m x ) = ^|^^<S^ SD (e, e). 

The convolutions of nuclear form-factors with the WIMP velocity distributions, A SI SD (e, e), are 
defined by expressions (22). We neglect for 73 Ge the sub-dominant contribution from WIMP- 
proton spin coupling proportional to (Sp 1 ). We consider only a simple spherically symmetric 
isothermal WIMP velocity distribution [31, 114] and do not go into details of any possible and 
in principle important uncertainties (and/or modulation effects) of the Galactic halo WIMP 
distribution [115-121]. For simplicity we use the gaussian scalar and spin nuclear form-factors 
from [113, 122]. For the relic neutralino mass density and for the escape neutralino velocity we 
use the values 0.3 GeV/cm 3 and 600 km/s, respectively. With formulas (25), we perform below 
a simple estimation of prospects for DM search and SUSY constraints with the high-spin 73 Ge 
detector HDMS assuming mixing of WIMP-neutron spin and WIMP-nucleon scalar couplings 
together with available results from the DAMA-Nal and LiXe experiments [1-3, 104, 111, 112]. 

The Heidelberg Dark Matter Search (HDMS) experiment uses a special configuration of 
two Ge detectors to efficiently reduce the background [9, 123]. From the first preliminary 
results of the HDMS experiment with inner HPGe crystal of enriched 73 Ge [9, 105] we can 
estimate the current background event rate R(e, e) integrated here from the "threshold" energy 
e = 15 keV to "maximal" energy e = 50 keV. We obtain _R(15,50) ~ 10 events/kg/day. A 
substantial improvement of the background (up to an order of magnitude) is further expected 
for the setup in the Gran Sasso Underground Laboratory. In Fig. 9 solid lines for the in- 
tegrated rate i?(15,50) marked with numbers 10, 1.0 and 0.1 (in events/kg/day) present for 
m wiMP = 70 GeV our exclusion curves expected from the HDMS setup with 73 Ge in the frame- 
work of mixed SD WIMP-neutron and SI WIMP-nucleon couplings. Unfortunately the current 
background index for HDMS is not yet optimized, and the relevant exclusion curve (marked 
with 10 events/kg/day) has almost the same strength to reduce <jg D as the dashed curve from 
the DAMA experiment with liquid Xe [23] obtained for m WIMP = 50 GeV (better sensitivity 
is expected with HDMS for m wlMP < 40 GeV). However, both experiments lead already to 
some sharper restriction for <7g D then obtained by DAMA (see Fig. 9). One order of magnitude 
improvement of the HDMS sensitivity (curve marked with 1.0) will supply us with the best 
exclusion curve for SD WIMP-neutron coupling, but this sensitivity is not yet enough to reach 
the calculated upper bound for crg D . This sensitivity also could reduce the upper bound for 
SI WIMP-proton coupling o"gj to a level of 10~ 5 pb. Nevertheless, only an additional about- 
one-order-of-magnitude HDMS sensitivity improvement is needed to obtain decisive constraints 
on (jgj as well as on (Xg D . In this case only quite narrow bounds for these cross sections will 
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FIG. 9: The solid lines (marked with numbers of i?(15,50) in events/kg/day) show the sensitivities 
of the HDMS setup with 73 Ge in the framework of mixed SD WIMP-neutron and SI WIMP-nucleon 
couplings. The DAMA-Nal region for sub-dominant SD WIMP-neutron coupling {0 = vr/2) is from 
Fig. 8. Scatter plots give correlations between cigj and <Tg D in the effMSSM for m x < 200 GeV. 
The squares (red) correspond to sub-dominant relic neutralino contribution 0.002 < fi x ^o < 0.1 an d 
triangles (black) correspond to WMAP relic neutralino density 0.094 < ^ x /iq < 0.129. The dashed 
line from [23] shows the DAMA-LiXe (1998) exclusion curve for rn^^p — 50 GeV. 

be allowed (below the curve marked by 0.1 and above the lower bound of DAMA-Nal mixed 
region). In practice it seems, that only the DAMA and the HDMS constraints together could 
restrict the SD WIMP-neutron coupling sufficiently. 



E. Some other consequences of the DAMA results 

It follows from Figs. 1, 5, 6 and 8 that the main results of the DAMA experiment one 
could summarize in the limitations of the WIMP mass, and the restrictions on the cross section 
of the scalar WIMP-proton interaction. Quite approximately (having in mind all possible 
uncertainties of [2, 3]) one can write them in the form: 

40 GeV < m WIMP < 150 GeV, 1 • 10~ 7 pb < cx^O) < 3 • 10" 5 pb. (26) 

The limitations of (26) should have some consequence for observables. Taking them into account 
we have obtained the reduction of our scatter plots for the total expected event rate of direct 
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FIG. 10: Event rate for direct neutralino detection in a Ge detector as function of the LSP neutralino 
mass. Green crosses present our calculations with relic density constraint 0.1 < ^ x /ig < 0.3 only. Open 
boxes correspond to implementation of the SI cross section limit MO" 7 pb < of r (0) < 3-10~ 5 pb only, 
and closed boxes show results with the additional WIMP-mass constraint 40 < w WIMP < 150 GeV 
(see (26)). 




FIG. 11: Indirect detection rate for upgoing muons from DM particles (neutralinos) annihilation 
in the Earth (a) and the Sun (b) as function of the LSP neutralino mass. Green crosses present 
our calculations with relic density constraint 0.1 < QyfiQ < 0.3 only. Open boxes correspond to 
implementation of the SI cross section limits of (26) only and closed boxes depict results with both 
limitations of (26). 

WIMP detection in a 73 Ge detector (Fig. 10) and the indirect detection rate for upgoing muons 
from dark matter particles annihilation in the Earth and the Sun (Fig. 11). The calculations 
of indirect detection rates follow the description given in [13, 124]. There is also a reduction 
of allowed masses of some SUSY particles (Fig. 12). In total from these figures one can see 
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m^, GeV m r GeV 

FIG. 12: Masses in GeV of light (m h ), heavy (?7%), and charged Higgs bosons (m H+ ), as well as 
masses of chargino (m + ), stop (?7i s top)> an d second neutralino (m 2 ), versus the mass (m x ) of the 
LSP neutralino under the same DAMA-inspired restrictions as in Figs. 10,11. 

that the DAMA evidence favors the light Higgs sector of the MSSM, relatively high event rate 
in Ge detectors, as well as relatively high up-going muon fluxes from the Earth and from the 
Sun for indirect detection of the relic neutralino. It is also almost insensitive to the sfermion 
and neutralino-chargino particle masses. As noted before in [44, 124] the relatively light Higgs 
masses (smaller than 200 GeV) are very interesting from the point of accelerator SUSY searches. 
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IV. CONCLUSION 

In the effective low-energy MSSM (effMSSM) for zero momentum transfer we calculated 
the LSP-proton(neutron) spin and scalar cross sections in the low LSP mass regime, which 
follows from the DAMA dark matter evidence. We compared the calculated cross sections 
with experimental exclusion curves and demonstrated that about a two-orders-of-magnitude 
improvement of the current DM experiment sensitivities is needed to reach the SUSY predictions 
for the <7go- 

We noted an in principle possible incorrectness in the direct comparison of exclusion curves 
for WIMP-proton(neutron) spin- dependent cross section obtained with and without non-zero 
WIMP-neutron (proton) spin-dependent contribution. On the other side, nuclear spin struc- 
ture calculations show that usually one, WIMP-proton (S^), or WIMP-neutron (S^), nuclear 
spin dominates and in the effMSSM we have the WIMP-proton and WIMP-neutron effective 
couplings a n and a p of the same order of magnitude (Fig. 7). Therefore at the current level 
of accuracy it looks reasonable to safely neglect sub-dominant WIMP-nucleon contributions 
analyzing the data from spin-non-zero targets. Furthermore the above-mentioned incorrect- 
ness concerns also the direct comparison of spin-dependent exclusion curves obtained with and 
without non-zero spin-independent contributions [2, 3]. To be consistent, for this comparison 
one has to use a mixed spin-scalar coupling approach (Figs. 8 and 9), as for the first time 
proposed by the DAMA collaboration [1-3]. We applied such spin-scalar coupling approach 
to estimate future prospects of the HDMS experiment with the neutron-odd group high-spin 
isotope 73 Ge. Although the odd-neutron nuclei 73 Ge, 129 Xe already with the present accuracy 
lead to some sharper restrictions for <7g D then obtained by DAMA, we found that the current 
accuracy of measurements with 73 Ge (as well as with 129 Xe and Nal) did not yet reach a level 
which allows us to obtain new decisive constraints on the SUSY parameters. Future about 
two-orders-of-magnitude improvement of the background index in the HDMS experiment [9] 
can in principle supply us with new constraints for the SUSY models. 

Finally we noticed that the DAMA evidence favors the light Higgs sector in the effMSSM 
(which could be reached at LHC), a high event rate in a 73 Ge detector and relatively high 
upgoing muon fluxes from relic neutralino annihilations in the Earth and the Sun. 

V.B. thanks the Max Planck Institut fur Kernphysik for the hospitality and RFBR (Grant 
02-02-04009) for support. 
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